A scalable, high-speed measurement-based quantum computer using trapped ions 
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We describe a scalable, high-speed, and robust architecture for measurement-based quantum- 
computing with trapped ions. Measurement-based architectures offer a way to speed-up operation 
of a quantum computer significantly by parallelizing the slow entangling operations and transferring 
the speed requirement to fast measurement of qubits. We show that a 3D cluster state suitable 
for fault-tolerant measurement-based quantum computing can be implemented on a 2D array of 
ion traps. We propose the projective measurement of ions via multi-photon photoionization for 
nanosecond operation and discuss the viability of such a scheme for Ca ions. 
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The rapid progress in quantum information processing 
systems has been fueled by the realization that the algo- 
rithmic complexity of a quantum computer scales polyno- 
mially with the size of certain important problems rather 
than exponentially [1 j. This presents a tremendous ad- 
vantage for large problems that are so far solvable only 
on a time scale of years. However, the practical util- 
ity of large-scale quantum computers will also depend on 
their ability to compete with current classical comput- 
ers. Consider, for example, Shor's factoring algorithm - 
the composite integer RSA640 (N = 640 bits) requires 
w 32 N 3 — 8.4 10 9 [2] operations (neglecting error cor- 
rection overheads and improved scaling with other trade- 
offs [3]). To compete with a distributed network that can 
factor RSA640 in 5 months [4], quantum operations on 
time scales of 1.5 ms are required; to factor RSA640 in, 
say, 5 minutes, operation time scales have to be improved 
to 36 ns. Error correction overheads will further worsen 
time scales. In other words, nanosecond operations are 
essential for practical large-scale quantum computers. 

In this Letter, we address this speed issue for one 
of the most promising quantum computing (QC) imple- 
mentations, namely, ion-traps [5 . In contrast to stan- 
dard quantum circuit schemes considered so far for ion- 
traps, we consider measurement-based QC paradigms [6], 
where the actual processor speed is mostly determined 
by the measurement time scales. We demonstrate (i) 
this one-way quantum computing (1WQC) scheme has 
significant advantages for ion-trap QC; (ii) that a 3D 
cluster state for fault-tolerant computing can be effi- 
ciently implemented in 2D ion-trap architectures; and 
(iii) that multi-photon ionization and detection of the 
emitted electron using the ion trap potential as a guide 
can significantly speed up high-efficiency measurements 
to nanosecond time scales. 

Ion- trap QC [5] has surpassed several major milestones 
on the QC roadmap [7]. Recent experiments have entan- 
gled up to eight ions [8] , demonstrated above 99% fidelity 
gates [9J[10], and coherence times of 10-34 s QTJ[12]. Dis- 
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tant entanglement of ions via interference of emitted pho- 
ton pairs has been demonstrated [13] and scalable chip- 
based trap architectures are being implemented p~4j [15] . 
However, time scales for logic-gate operations are slow, 
on the order of 1-100/is for entangling gates [5j [16] and 
1/as to 10ms for single-qubit operations and measure- 
ment [10 . Moreover, the shuttling of ions required in 
scalable architectures puts even worse timing constraints 
on two-qubit gates (50-100/is) [17 . Measurement-based 
QC paradigms [6] offer a way around this because the 
computational resource, a multipartite entangled cluster 
state, can be created via entangling operation applied in 
parallel and offline. This has the tremendous advantage 
that the usual requirement to avoid motional heating in 
ion traps is restricted to this first entangling step and is 
removed from the actual QC process. In the 1WQC pro- 
tocol, processing of information occurs through measur- 
ing qubits in a prescribed basis, combined with feedfor- 
ward of measurement outcomes. The processor speed is 
determined by measurement, readout, and classical pro- 
cessing time scales which, as we show, should be possible 
on a nanosecond time scale. Methods for error correction 
have recently been proposed which introduce the neces- 
sity for 3D cluster states [18] . We describe how a 3D 
cluster can be efficiently implemented in a 2D architec- 
ture without resorting to time as third dimension [19] . 

In the 1WQC paradigm, all entanglement operations 
are done in parallel and offline before commencement 
of the algorithm. This multipartite entangled cluster 
state is created by shuttling ions in a 2D lattice ge- 
ometry into close proximity and by applying standard 
controlled-PHASE (CPHASE) gates [3 US [20] between 
neighboring atoms. In such a 2D array, measurements of 
the qubits in different bases and feedforward of measure- 
ment outcomes allow the simulation of a universal quan- 
tum circuit [6]. While this 2D cluster state represents 
a universal resource, topological encodings in 3D struc- 
tures are required for error correction with high thresh- 
olds [18] . More recently, quasi-2D fault tolerant cluster 
states, where the third dimension is replaced by time, 
were introduced [19] , effectively mapping a 3D cluster to 
a 2D array. However, in this scheme, the two-qubit gates 
required during the measurement-based QC will nullify 
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the potential for speed-up over circuit based QC. 

Because of the requirements of laser access, implement- 
ing a full 3D architecture for fault- tolerant 1WQC using 
ion traps would be extremely challenging. Instead, a 3D 
cluster state can be implemented efficiently in a 2D ion- 
trap array with non-nearest neighbor entanglement op- 
erations. The architecture required for implementing the 
1WQC paradigm is similar to the scalable ion- trap archi- 
tectures proposed in [14 . This architecture requires in- 
herently slow shuttling of ions [17 . For faster transport, 
better fidelity, and to keep ions in the motional ground 
state, 120° Y-junctions in a 2D hexagonal array (Fig{TJ) 
are preferable to 90° junctions in a square lattice [T5] . 




FIG. 1: (a) Hexagonal ion trap architecture with 120° junc- 
tions. Single channel electron multiplier (CEM) detectors for 
high-efficiency electron detection are indicated. For entan- 
gling operations, atoms can be shuttled into close proximity 
with high fidelity to remain in the ground state. To entangle 
non- nearest neighbor atoms, one atom could be moved out 
of the way as indicated by arrows. Alternatively, swap op- 
erations between qubits at a junction could minimize move- 
ment of atoms, (b) Underlying hexagonal array and result- 
ing cluster states for non-nearest neighbor interactions. The 
hexagonal array is decomposed into two rhombic lattices (red 
solid and blue dashed lines show the edges of the cluster state 
layers). Each lattice can again be broken down into further 
sublattices. For example, 4 sublattices of the solid line rhom- 
bic lattice are indicated by circles, squares, triangles and dia- 
mond. These, together with 4 sublattices of the dashed rhom- 
bic lattice can create an eight-layer 3D cluster state (see text 
for details) . Qubits belonging to the eight different layers are 
numbered from 1 to 8. (Color online.) 

A full 3D cluster state can be created in this univer- 
sal [21] hexagonal 2D architecture in the following way: 
If we do not restrict ourselves to nearest neighbor entan- 
gling operation, then the hexagonal lattice can be broken 
down into two rhombic lattice sublattices with distance 
2d between sites as illustrated in Fig[]Jb). Each of these 
two sublattices then represents a single 2D layer of a 
3D cluster state. To create each layer of cluster state, 
atoms in each sublattice have to be entangled with its 
sublattice neighbors, 2d away, as described in Figjlja). 
In each sublattice, a sequence of four CPHASE gates ap- 
plied in parallel across the array entangles every ion with 
all four sublattice neighbors, creating a full 2D cluster 
state layer [6] . Entanglement between different layers can 



be accomplished by simply entangling ions belonging to 
different sublattices (indicated by numbers 1-8) in just 
two parallel applications of a series of CPHASE gates as 
indicated by thick black lines in Fig. 1(b)]. The number 
of layers can be increased from 2 to 2n 2 by increasing 
the elementary cell of each sublattice by a factor n. This 
corresponds to entangling atoms separated by a distance 
2nd [see Fig. 1(a)]. Note that the distance only scales as 
2n whereas the number of layers is 2n 2 . An example of 
an 8-layer system is shown in Fig. 1(b) for a distance of 
4d between entangled atoms. Another important feature 
is the possibility to entangle the first and last layer of 
the 3D cluster state, thereby directly and efficiently cre- 
ating a topological structure as required in [THJ [19] . The 
number of parallel operation and thus the time to create 
this 3D cluster state offline is constant in the number of 
qubits, and only requires six CPHASE gates cycles ap- 
plied in parallel to all neighboring qubits. 

In addition to standard entangling operation between 
nearby ions, probabilistic distant entangling operations 
[13] [22] could be used to connect 3D "sub" -clusters, 
which are created in separate locations and contain 10 3 - 
10 4 qubits. Building up a large cluster this way 1) avoids 
correlation errors across the entire cluster state, 2) pro- 
tects the unmeasured cluster from disturbance due to 
the ionization measurements, and 3) allows for slightly 
delayed preparation to avoid storage errors as well as 
re-preparation of subclusters while others are measured. 
This fulfills the requirements of the more stringent error 
model in [18 , which includes errors due to correlations 
and storage (error model 2 [H]). Due to disparate time 
scales between measurement (~ns) and decoherence time 
(~s) storage errors are kept well below 10 -4 even when 
measuring up to 10 4 qubits. In addition, decoherence 
free subspaces have been successfully demonstrated for 
ions [12] as well as for 1WQC [23] . A detailed error anal- 
ysis with accurate thresholds and overheads is planned 
for future work. 

Measurements of atomic qubits are usually achieved in 
two steps: A single-qubit rotation, followed by projection 
into one of the two-qubit states. We assume encoding 
in the optical S-D transition in Alkali- metal- like ions to 
facilitate nanosecond single qubit rotations and readout 
using light pulses with a bandwidth much smaller than 
the energy separation between the S and D states. To 
increase decoherence time a mixture of encodings would 
be conceivable, i.e. long-time storage in the hyperfine 
state and short-time storage in S-D states for faster ma- 
nipulation. Transitions between the S ground state and 
metastable D state are electric-dipole forbidden with long 
coherence times (~s). Single-qubit rotations on the S-D 
transition can be achieved via quadrupole-allowed tran- 
sitions or two-photon transitions using the intermediate 
P\/2 or P3/2 states. Quadrupole Rabi frequencies Slf^ 
in current experiments are 35.5 kHz for fairly low irra- 
diance of about 6W/cm 2 [10 . To achieve manipulations 
on a 2 ns time scale via quadrupole fields would require 
extremely high ir radiances of 10 9 W/cm 2 . For Raman 
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transitions, detuned 10 4 line widths from the P states, 
an irradiance of 10 5 W/cm 2 is sufficient to achieve ma- 
nipulations on a 1 ns time scale. 

A more severe timing limit and one that requires a 
more drastic shift in approach, is due to the fluores- 
cence readout in ion systems where one of the two qubit 
states is coupled to a shelving state via a cycling tran- 
sition [10]. The time scales for this are limited by the 
lifetimes of dipole transitions, and photon-collection and 
photodetection efficiencies to a few microseconds. How- 
ever, state-dependent multiphoton ionization and sub- 
sequent detection of the electron allow measurement on 
nanosecond time scales. Consider a resonant four-photon 
ionization of Ca + in the 5-state (see Fig. |2| using easily 
accessible 380 — 410nm excitation wavelengths. From the 
level diagram in Fig. [2j it is immediately apparent that a 
multiphoton ionization transition from the 4,Si/ 2 -ground 
state should be possible for fairly low intensities. On 
the other hand photoionization of the D-level is strongly 
suppressed, as the detunings are an order of magnitude 
larger for each transition. For 380 — 410nm, the four- 
photon transition is close to resonance with several levels 
for ionizing atoms in the 4S , 1 / 2 state. We can choose 
a resonance condition with the 4P x / 2 state at 397nm, 
the 551/2 state at 383nm, or the 6Pi/ 2 and 6P 3 / 2 states 
at 403nm. A broadband frequency-doubled Ti: Sapphire 
with appropriate pulse shaping would be able to address 
all three resonances and exploit interference between ion- 
ization paths to improve the ionization fidelity and state 
discrimination l24l. 



only near-resonance transitions in our estimate of the ion- 
ization rates [26^ (non-listed matrix elements have been 
estimated by simple scaling). Using this basic estimate, 
we expect the nonresonant parts of photoionization prob- 
ability to be a factor of at least [( J*)s/( J*)d] 2 = 1600 
lower for atoms in the D state compared to the S state. 
This difference allows a very high fidelity state discrim- 
ination as would be necessary for a projective measure- 
ment. At resonance with the 6Pi/ 2 state, we expect ion- 
ization rates on the order of « 10 9 — 10 10 s _1 for 
peak irradiances of about / = 10 9 W/cm 2 . This allows 
for ionization on a femto- or picosecond time scale with 
accessible laser intensities 24 . 
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FIG. 2: Relevant levels for ionization-readout scheme of the 
(a) 5-state and (b) D-state in 43 Ca + . Contributing transi- 
tions are shown with black arrows. The ionization threshold 
at 11.87 eV is indicated as thick solid line. (Color online.) 

Using the effective operator method [25 , the reso- 
nant ionization probability for the 5-state when ap- 
plying laser light of irradiance / and tt polarization is 
N% = E A 47r/4 ( J 7?) 2 +47r/ 2 K 2 /L 2 . The nonresonant 
contributions are described by transition probabilities 
for transitions to continua states with different angular 
momentum A = 5, P, D, etc.; L, K are resonant tran- 
sition operators [25]. For the nonresonant D-state ion- 
ization, only an incoherent sum of the transition rates, 
= 47r(J^ ) 2 / 4 need be considered. We include 



FIG. 3: (a-c) Time evolution of the emitted electron 
wavepacket in the ion-trap saddle potential at three different 
times, (a) 0.25 ns (b) 1.0 ns (c) 1.5 ns. We assume an initial 4- 
photon kick of the electrons in the positive x-direction. (d) If 
we place one detector 30/xm away on the x-axis, we can detect 
the electron with up to 90% fidelity [see dashed line in (d)]. 
A second detector can be placed on the negative x-axis to de- 
tect the part of the wave packet that escaped along the saddle 
in this direction. The detection efficiency is then above 99% 
within a 3 ns time window [solid line in (d)]. (Color online.) 

Photoionization and detection of ions has been used 
in multiple experiments, e.g. for counting of atoms [27] , 
imaging of a gas of atoms [28] and even single atoms in 
a cloud [29]. Detection of electrons via a single chan- 
nel electron multiplier (CEM) is generally more efficient 
(> 99% efficiency) than detection of the heavy ion [30] . 
A CEM also features detector dead times of nanoseconds 
or less [30]. Not only is the detection probability higher 
for the freed electron compared to the ion, but electrons 
move faster, so that the time between ionization and de- 
tection is less. Stability of a particle in the ion trap 
requires usually the characteristic trap parameter q x to 
be in the range < q x < 1. The electron is unsta- 
ble in the dynamic trapping potential as q e x oc ^e/m e 
is two orders of magnitude larger than the one for the 
ion <7^ on oc ^/ e / m ion- The rf frequency ou v { for the ion- 
trap field is typically on the order of 20-30 MHz [10] so 
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that the characteristic time scale for the electron motion 
to « l/^ r fy/m e /mion in the ion trap is about 0.5 ns. 
The ion-trap field is therefore time-independent on the 
electron detection time scale, i.e., the electron moves in 
an almost time-independent saddle-point potential and 
will be ejected along the saddle. To describe the mo- 
tion of the electron in this potential, assume all momen- 
tum from the four photon ionization is transferred to 
the electron in the positive x-direction (the inertia of the 
trapped ion is effectively infinite); the corresponding ini- 
tial electron velocity is then 7 x 10 3 m/s. The electron is 
ejected when the dynamic trap potential is anti-trapping 
in the x-direction and trapping in the y-direction (with a 
static weak trap field in the longitudinal z-direction [TO]). 
A time-dependent simulation of the electron wavepacket 
dynamics shows that the ion-trap potential can act as 
a guide for the electron onto the CEM detectors (see 
Fig. 3). Our simulations show that the electron remains 
localized in the y-direction with the typical wave packet 
breathing due to squeezing in the y-direction. In the x- 
direction, the electron wave packet spreads and "slides" 
down the saddle of the trapping potential due to the ini- 
tial momentum kick of the photoionization photons. If 
we place two 20/im size CEM detector about 30/im away 
from the trap center close to the trap-electrodes on the 
positive and negative x-axis, we can estimate above 99% 
detection efficiency within 3 ns [Figj3^d)]. Due to trap 
geometry, placement of detectors may be restricted, but 
the initial conditions (direction of photoionization beam, 
timing of ionization) should make it possible to guide the 
electron into the detector using the existing trap-field. 

For above ionization readout, we also need to con- 
sider what happens to the remaining doubly ionized Ca 2+ 
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meters apart [31]. For quantum computing, a detailed ex- 
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accurate error correction thresholds and overheads. 
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